Updated constraints on f(7l) from cosmography 
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We address the issue of constraining the class of f(JZ) able to reproduce the observed cosmological 
acceleration, by using the so called cosmography of the universe. We consider a model independent 
procedure to build up a /(2:)-series in terms of the measurable cosmographic coefficients; we therefore 
derive cosmological late time bounds on f{z) and its derivatives up to the fourth order, by fitting 
the luminosity distance directly in terms of such coefficients. We perform a Monte Carlo analysis, 
by using three different statistical sets of cosmographic coefficients, in which the only assumptions 
are the validity of the cosmological principle and that the class of f{TZ) reduces to ACDM when 
0^1. We use the updated union 2.1 for supernovae la, the constrain on the Hq value imposed by 
the measurements of the Hubble space telescope and the Hubble dataset, with measures of H at 
different z. We find a statistical good agreement of the f{TV) class under exam, with the cosmological 
data; we thus propose a candidate of f{TZ), which is able to pass our cosmological test, reproducing 
the late time acceleration in agreement with observations. 

PACS numbers: 04.50.+h, 04.20.Ex, 04.20.Cv, 98.80.Jr 



I. INTRODUCTION 

The recent observational evidence of the late time ac- 
celeration of the universe [l|, [1| opened new challenges 
in the framework of theoretical cosmology. To explain 
the origin of such a cosmic speed up, cosmologists usu- 
ally assume the existence of an exotic fluid called dark 
energy (DE) Q. Even though its physical nature is still 
unclear, several attempts have been made in order to re- 
solve the problem of its existence Q . In general relativity 
(GR), the simplest possibility is the introduction of a vac- 
uum energy cosmological constant, A [^,0]. The resulting 
model is usually referred to as ACDM Q. However, al- 
ternative approaches have followed each other, without 
being conclusive (8l-[ll|. To this regard, another appeal- 
ing possibility is to consider GR as a limiting theory of a 
more general paradigm [l^]; so that, in the last decades, 
particular attention has been devoted to solve the prob- 
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lem of the accelerated universe in the framework of ex- 
tended theories of gravity [l3|. Generally, extending GR 
means to review the DE effects as due to possible correc- 
tions of the Einstein- Hilbert action • 

In this paper, we focus our attention to the case of the 
so called f{Tl) theories, in which the Ricci scalar TZ in 
the Einstein-Hilbert action is replaced by a more general 
analytic function, namely f{TV}. The corresponding ac- 
tion reads A — j d^x^/^ [/('^) + -Cm] [II], where Cm is 
the standard matter term. By varying the action with 
respect the metric , we obtain the field equations [T6j 

T^M.^/' W-^/(7^)ffM'^-(V^V,-5^,VaV")/'(7^) = SttT^ 

(1) 

where T^^ is the standard energy momentum tensor and 
G = c=l. 

The problem of determining the nature of DE is there- 
fore shifted to understand which f{TV) is the correct can- 
didate to explain the dynamics of the universe. The cor- 
rect class of fjJZ) should be compatible with modern 
observations [13] ■ Therefore, we propose to limit our at- 
tention only to the class oi f(TZ) reducing to ACDM at 
the low redshift regime |18l - [23 1. 

In order to check the viability of f{Tl) models, it has 
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been proposed in {23j to study the so called cosmography 
of fijV). Cosmography represents a part of cosmology 
which does not postulate any cosmological model a pri- 
ori. Thus, it can be thought as a model independent way 
to fix constraints on the universe dynamics at late times 
through the use of a set of parameters, namely the cos- 
mographic set (CS). The aim of this work is to relate the 
f{TV} Taylor expansion around z ^ to the luminosity 
distance, determining the Taylor coefhcients as functions 
of the CS. Afterwards, we fit the coefficients by directly 
using the luminosity distance; this allows us to overcome 
the problem of the error propagation, since the /(??,) co- 
efficients are measured directly from data. In particular, 
once obtained the expression of the luminosity distance 
in terms of f{z) and its derivatives, we perform a Monte 
Carlo fitting procedure. We obtain at our time strin- 
gent numerical intervals for f{z) and its derivatives up to 
the fourth order and then the corresponding constraints 
on f{TZ) and its derivatives. The set of cosmographic 
f{z) parameters is measured directly from supernovae la 
(SNela) data, 'H(z) observations, under the bound pro- 
vided by the Hubble space telescope for Hq. Finally, we 
get a viable candidate of f{Tl), reconstructing it from 
the cosmographic test. Such a candidate seems to pass 
the cosmographic and cosmological tests, extending the 
ACDM paradigm as a limiting case of a more general 
theory. 

The paper is structured as follows: in Sec. II we de- 
velop the main features of cosmography and we define 
the so called f{z) cosmographic set, which is the set of 
f{z) and its derivatives to be fitted with the data. In Sec. 
Ill we perform a Monte Carlo analysis, based on three 
statistical models, while in Sec. IV we propose a viable 
candidate of f{Tl), compatible with the bounds inferred 
from our tests. Finally, in Sec. V we develop conclusions 
and perspectives of our work. 



II. COSMOGRAPHIC f{z) PARAMETERS 



[26|. We give here the expressions for the first 6 coeffi- 
cients in the expansion. 



I = 



a dt 
1 d^a 
" aU^ dt^ 
I d^a 



1 



T,W> dt5 



aH^ dt^ 
1 d^a 
aU^ dt^ ' 
1 d^g 



(3) 



The coefficients in eq. ^ are, by construction, model in- 
dependent quantities, which are called the cosmographic 
set (CS). They are known in the literature as the Hubble 
rate (H), the acceleration parameter (q), the jerk param- 
eter (j), the snap parameter (s), the lerk parameter {I) 
[m and the m parameter introduced in (22]. The set of 
such parameters is known in the literature as the CS. 



A. Degeneracy and cosmography 

The definitions given above lead to the most rele- 
vant property of cosmography, that is, overcoming the 
so called degeneracy problem among different cosmolog- 
ical models. In fact, no cosmological model is assumed a 
priori in the expression of the luminosity distance. Fur- 
thermore, another significative aspect of cosmography is 
to relate the series expansion of the luminosity distance 
to the CS. To this regard, it was pointed out |24l - [26| that 
direct measurements of such quantities are permitted, 
overcoming the problem of the statistical error propaga- 
tions. Hence, it is possible to compare theoretical predic- 
tions with the observed values, without passing through 
a cyclic scheme which postulates a priori the form of T-L 
and /(i?) [H. 

One of the most important observational quantities to 
be expanded in series is the luminosity distance d^. Con- 
sidering the scale factor definition in terms of z, i. e. 
a = (I + z)^^, the luminosity distance reads 



In this section, we relate the f{TV) coefhcients (evalu- 
ated in terms of the redshift z) to the cosmographic set 
(CS). Afterwards, we use these relations to write the lu- 
minosity distance in terms of /(z) and its derivatives at 
z = 0. To this end, let us review briefly the theoretical 
features of cosmography. Cosmography, or alternatively 
cosmo-kinetics, is a tool to investigate the dynamics of 
the universe, regardless of the particular cosmological 
model. Cosmography indeed simply postulates the va- 
lidity of the cosmological principle. Thus, it follows the 
use of the Friedmann-Robertson- Walker (FRW) metric, 
i.e. 



ds^ 



dt^ - a{tf {dr^ + r^dn^) 



(2) 



where we assume hereafter a spatially flat universe (fc — 
0) and we use the notation dil.^ = d0^ + sin^ OdcjP . 

The paradigm of cosmography is to expand the scale 
factor a{t) in a Taylor series around the present time to 



dL = 



C 

47rJ' 



_ro_ 

ait)' 



(4) 



where we defined C and T as the luminosity and the flux 
respectively, while 



_d^ 

«(0 



(5) 



whose physical meaning is related to the distance r that 
a photon travels from a light source at r — ro to our 
position at r = 0. Equation (j4|) can be expanded in 
powers of z around z — 0; the expansion up to the sixth 
order in z is given in the Appendix A, both in terms 
of the CS and in terms of the derivatives of /(z). Now 
we want to write f{TZ) — f{JZ{z)) = f{z) and use the 
definitions in 1^ to express /(z) in terms of the CS, i.e. 
/(z) = /(H(z),q(z),j(z),s(z),Z(z),m(z)). 
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To do so, let us start from the definition of TZ in terms 
of t and i.e. 

n = -G(n + 2H^^ . (6) 
Using tlie redshift definition in terms of the cosmic time 



dlog(l + z) 
dt 



we rewrite TZ in terms of z, obtaining 

n = 6[{l + z)HH,~2H^] 



(7) 



(8) 



Hence, we can calculate TZ and its derivatives in terms of 
z and evaluate them in z = 0. The result, up to the fifth 
derivative, is 



TZo 
6 

TZzo 

6 
6 

6 

Ti-5z0 

6 



: [HzO ~ 2Ho] I 

: — 2^20 + STizo'H-2zO + Tio{ — 2Ti2zG + Tiszo) : 

' 3^220 + ^zo (~3'H2zo + 4^320) + 'Hoi—Tiszo 

+ ^4zo) J 

: WTi2zoT~l-3zQ + ^'Hz^yHizQ + 'HiyH.^zO , 
: 10'H32o('H220 + T~izza) + 15H2z0^4zO 

+ 'Hzo (5'H4zo + 6^520) + Ho {Tibzo + Hezo) , 

(9) 



where, hereafter, we adopt the convention 4^ = ^nzOj 

for X a generic function of z. 

Therefore, in order to evaluate TZ = 
TZ{T-Lo,qo,jo,so,lo,mo), we need to express Ti and 
its derivatives in terms of the CS. To this regard, after 
some cumbersome algebra, we infer from Eqs. ([3]) 



1 = -^-1 



^ 1/3 



3q-2, 



^ + 4j + 3q + 4) + 6 , 

7^(4) 



(10) 



7^(5) 



- 24 - 60(7 - 30g2 - lOj (g + 2) + 5s , 
+ 10j2 + I20j (<7 + 1) + 



3 [2/ + 5 (249 + 18g2 + 2g^ - 2s - gs + 8)] , 

and then the corresponding derivatives of Ti in terms of 
the cosmic time read 



=Ti^ [s-4j-3g(<7 + 4)-6] , 

[I - 5s + lQ{q + 2)j + 30{q + 2)q + 24] 



d^n 

dt^ 

=-H^|m - 4? + 12s + Isq - 24j - 32jg - Af 

- 2Aq - 36g2 - 6q^ - 6(j + 3q + 2f 
+ 8(l + g)(s-4j-3g(g + 4)-6) 

- 2[; - 5s + (lOj + 30g)(g + 2) + 24] } . 

(11) 

Thus, using Eq. ([7]), we can rewrite Eqs. ([TT|) in terms 
of the CS only, obtaining 

T~Lzo/Tin = 1 + 90 I 
'H2zo/'Ho =io — 90 J 

'H^zq/'Ho = - 3jo - 4^090 + 9o + S^o " ■'o , 

H4Z0/H0 = 12io - 4j2 + ^0 + 32jogo - 12(Z^ + 25^090 

- 24^3 - I5g,4 + + 7^p5g ^ 

Hszo/Ho = 32jogo - 15^0 - U^ogo + 60g^ + 180g^ 
+ 225g^ + 105g^ + 10io'(6 + 7go) - mo 

- jo(60 + 272^0 + 375^2 + 2lQql - 15so) 

- 60so - 98goSo - 60goSo - 7goSo • 

(12) 

Then, using equations Q and we are able to eval- 
uate the expressions of TZ and its derivatives as functions 
of the CS only. 

B. The use of the modified Friedmann equations 

In this subsection, we want to show the procedure to fix 
constraints on f{TZ) and its derivatives. We therefore use 
Eqs. (O and IT^ and we consider the modified Fried- 
mann equations, derived by assuming the FRW metric 
and Eq. Q. 

In the case of the standard matter term, (p,„ cx 
and P,n = 0) , one gets the modified Friedmann equations 



Ti' 



and 



-\- 31~L — Pcurv ■ 



(13) 



(14) 
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Equations p3)) and ([14]) determine the definition of tlie 
DE fluid in terms of the curvature as 



1 f 1 



/'(7^) i2 



fin) - 7^/'(7^) 



- mnf'in) 



The corresponding barotropic pressure reads 

Pcurv — ^curv Pcurv ; 



(15) 



(16) 



with the definition of the effective curvature barotropic 
factor given by 



= -1- 



iif" (TZ) + n nf"'{n) - uf'in) 



[/(7^) - 7^/'(7^)] /2 - snnfin) 



(17) 



Assuming the functional dependence TZ = Tl{z), we 
rewrite each term of Eq. pT)) in terms of z. We get 



/'(7^)=7^;l/., 
/"(7^)-(/2.7^,-^7^2.)7^,-^ 



/"'(7^) - 



(18) 



7^5 ' 



and, using equation ([7|), 
7t = - (1 + z)-H7e3 , 

a =(1 + z)H [H7^. + (1 + z){'Hz'Rz + H7e2.)] . 



(19) 



Furthermore, fohowing [23|, we know that any fiJZ) 
theory requires 



/"(7^o)=0, (20) 
in order to be compatible with Solar System tests and 



/'(7^o) 



1, 



(21) 



to predict the correct value for the gravitational constant 
G. 

Therefore, combining equation with ^ and (IT^ . 
we have 



ml 

fzO 

ml 

f2z0 _ TZ2 

ml ~ 



nr. 



qo 



TZzo _ „ 



-Jo 



(22) 



GUI 



-2 - 4qo - (2 + qo)jo ~ So : 



/(Ko)-Ko 

2 



and 



where we used the condition that Pcurvo 
that /o = 6H§(1- f7„) +7^o. 

Now, using equations (IT51) and ((T^ in p7|) . we can 
write Wcuri; as a function of z only. Then we expand this 
expression as a Taylor series around ^ = 0, obtaining 



E 

J=0 



1 d^c 



2 = 



(23) 



The first term in this expansion, which we call wo, 
depends only on / and its derivatives up to the third 
order (evaluated at ^ = 0), while the second term ioi 
depends on / and its derivatives up to the fourth order, 
and so forth for the higher terms. 

Keeping in mind that the class of f{Tl) should reduce 
to ACDM at low redshift regime, we compare our results 
with ACDM; thus, by fixing in equation ^ the ACDM 
bounds 



= 0(1 + 90) 



ACDM 
'^O 



,ACDM 



-3(1-290): 



(24) 



we get /o, fzo, hzo, fszo and /4^o in terms of the CS 
only 



fo 

2Hl 

f: 



zO 



ml 

f2z0 

ml 

2Hl 



J4z0 

2Hl 



qa 



= - 2 - qa + jo , 



- 2 - 4(70 - (2 + qo)jo - sq , 

- 4 - 3j2 + 3lo + Jo [2 + (7o(13 + 5qo)] 
f 15so + qo[2 + 2qo{h + 2qo)+ 9so] , 

+ 30jo(l + qo) - 6/0(5 + iqo) - 3mo - 66so 

- jo [22 + go (46 + go (38 + 29go)) - 15so] 

- go [18 + 84so + go(4 + 2qo{-Q + 2qo) 
+ 33so)] . 

(25) 

We refer to Eqs. as the definition of the f{z)- 

cosmographic set (fCS). Now our intent is to constrain 
the values of fo, fzO, /2z0, /32O and f4zO- To do so, we 
write the luminosity distance in terms of the fCS by 
using Eqs. (j25p . This is performed in two steps; first we 
invert the algebraic system ((25|) to find the CS in terms 
of the fCS. Then we insert these expressions in equation 
(|A1[) . The result is as a power series of z, whose co- 
efficients are now in terms of the fCS, instead of the CS. 
The explicit expression of in terms of the fCS is given 
in Eq. ([M]) . 

In addition, in order to measure the fCS using cIl and 
the cosmological data, we need to define viable priors. 
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Flat priors 



0.5 


<; 


h 




0.9 


n nni 








n riQ 


-5 


< 


W-\fo 


< 


5 


-10 


< 


10~V.o 


< 


10 


-15 


< 


IO-V2.0 


< 


15 


-20 


< 


10-^3.0 


< 


20 


-50 


< 


10-^4.0 


< 


50 



Additional constraints 

^^fc = 
= 0.274 
= 



Table I: Priors imposed on the parameters in the Monte Carlo 
analysis. 



compatible with the observed universe. To infer our pri- 
ors we assume that the class of /(??.) reduces to ACDM 
at late times, as already stressed above. We write such 
priors in Tab. I. 

We can now perform a best fit for the values of the 
fCS and obtain constraints on the values of /(z) and its 
derivatives at present time. This will be the content of 
the following section. 



• The union 2.1 SNela compilation of the supernova 
cosmology project [271]. This sample is an update 
dataset of the previous compilations union 2 [28j 
and miion 1 29] . Union 2.1 includes measurements 
in the plane /i— z of 580 supernovae over the redshift 
range 0.015 < z < 1.414. In the following numeri- 
cal analyses, we take into account systematic errors 
in the covariance matrix. 

• Observations of the Hubble factor (OHD) as a func- 
tion of redshift. We take the compilation of refer- 
ence [soj which encompasses 18 measurements be- 
tween the redshift range 0.09 < 2 < 1.75 (see Tab. 
I of [l^l). The data are extracted from previous 
works (see for ex. [3ll - l33| ). 

• A gaussian prior on the Hubble constant of Ho = 
74.2±3.6km/s/Mpc [H], as measured by the Hub- 
ble Space Telescope (HST). 

To constrain the parameters, we use a Bayesian 
method in which the best fits of the parameters are those 
which maximize the likelihood function 

£ (X exp(-xV2) , (29) 

where is the {pseudo)chi-squared function [ssj . 
Since the different sets of observations are not correlated, 
the function is simply given by the sum 



III. MONTE CARLO ANALYSIS AND 
CONSTRAINTS ON FCS 

In this section we evaluate the cosmological constraints 
on the fCS by fitting the luminosity distance (IA2p with 
the cosmological data. We analyze three statistical mod- 
els with different maximum order of parameters; this pro- 
cedure, widely adopted in the literature, corresponds to 
assume a hierarchy among parameters. The sets that we 
are going to analyze are summarized as 

A = {Ho,./o,/.o,/2.o}, (26) 
B = {Ho, fo, fzo, f2zo, fszo} , (27) 

C = {Ho, fo, fzO, f2z0, fszO, f4zo} ■ (28) 

In particular, the reason for studying the fCS in such 
a hierarchical way is that it is naively expected a broad- 
ening of the sampled distributions by adding more pa- 
rameters. The corresponding numerical effects to the 
measured quantities lead to strong error propagations; 
this is due to the higher orders of the Taylor expansion. 
We are interested both in quantifying these effects and 
in fixing constraints on the fCS. Our numerical study is 
based on a Monte Carlo simulation, in which the chosen 
observational datasets for our fits can be summarized as 
follows 



= Xu„ion2.1 + XhST + XOHD ■ (30) 

We perform a Markov Chain Monte Carlo analysis by 
modifying the publicly available code CosmoMC [3Q| (see 
also [S^l). To obtain the posterior distributions, we as- 
sume uniform priors over the intervals given in Tab. I. In 
Tab. [U we show the summary of the constraints. We re- 
port the best fits given by the maximum of the likelihood 
function of the samples, the quoted errors show the 68% 
confidence level (c.L). In Fig. 4 we plot the correspond- 
ing posterior distributions. The vertical lines denote the 
upper and lower limits for the ACDM case, these are 
obtained by using the best fits parameters reported in 
Tab. I, compatible with those of [s^. In Figs. 1, 2 and 
3, we show all the 2-dimensional marginalized posterior 
confidence intervals for the three analyzed models. 

As it can be noticed from figures 1, 2, 3 and 4, the 
marginalized posteriors loose Gaussianity when we add 
further parameters to Model A. We conclude that con- 
sidering Model C over Model B has the advantage that 
it gives more information on the cosmographic f(TV) 
parameters without enlarge the dispersions; however. 
Model C is less suitable for a posterior statistical treat- 
ment. 

We note that the dispersions of the samples are con- 
siderably enlarged when the third derivative fszo is in- 
cluded within Model A. In other words, the correspond- 
ing Model B suffers from a deep dispersion problem due 
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Table II: Best fits of the parameters for the tfiree considered models. Tiie quoted errors show the 0.68 c.l. The observations used 
to constrain the parameters are the union2.1 data set compilation, Observational determination of the Hubble factor (OHD), and the 
measured value of Ho by the HST. 



Parameter 


Model A 
xLn = 529.0 


Model B 

Xl^^r, = 540.0 


Model C 

Xmin ~ 552.6 




77 9*?+o.S4 


yc f;q+2.03 


71 'Hn+i 92 


10-Vo 


-3.324±S-- 


— 3 14.4+0-320 

^•-^^^-0.332 


-2.669iS-- 


10-V.o 








lO-Va.o 




2-27611111 


n 400+2.192 

u.-±aa 2.049 


10-^3.0 




-8.264+1-™/, 


—0 399+1 424 


10-^4.0 






— 1 027+2-43" 
J^-'-'^' -3.132 



Notes. "Ho is given in Km/s/Mpc. 



Table III: Standard proportions of f{z) derivatives. 



Parameter 


Standard deviations proportions 




Model A : Model B : Model C 


■Ho 


1 : 1.51 : 1.44 


fo 


1 : 1.43 : 1.25 


/zO 


1 : 3.26 : 2.59 


f2z0 


1 : 2.26 : 2.19 



Figure 1: 2-dimensional marginalized probability for the pa- 
rameters of model A. The dashing denotes the likelihood of 
the samples. 



Figure 2: 2-dimensional marginalized probability for the pa- 
rameters of model B. The dashing denotes the likelihood of 
The samples. 



IV. EXAMPLES OF /(7^) GRAVITY 



to the considered dataset of 580 SNela. Nonetheless, the 
introduction of /420 in Model C does not substantially 
broaden the distributions. To quantify such an effect, 
the standard deviations of the distributions are in the 
proportions given in Tab. IIIII An additional comment 
comes from the strong tension between Model A and the 
ACDM model; such a tension could be substantially al- 
leviated by considering Model B and Model C. 



In this section, we provide a new explicit example of an 
f{Tl) model that reduces to ACDM when 2: ~ and sat- 
isfies the theoretical constraints (|20|) and (f2T|) . Further- 
more, we set the free parameters of the model according 
to the new constraints on higher order derivatives that 
we found in Sec. Ill from cosmography. The expression 
for our f{TZ) candidate is derived in accordance with the 
above results, through the inverse procedure of determin- 
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Table IV: Values of the cosmographic set. 



Parameter 


Model A 


Moael J3 


Moael u 


go 


u. ( c^^_o,324 


— n 744+" 'i26 

w. 1 '±'±_o.434 




jo 




W.O-L 1 „2.102 


0.787tlit 


so 






-2.217tll°3 


lo 








mo 






-41.781iti;f3 



Table of numerical results for the CS; the numerical values are given at ^ = 0, while the error propagations have been found using the 

standard differential rule. 



Table V: Values of f{R) and its derivatives. 



Parameter 


Model A 


Model B 


Model C 




— S 'iOA+0.227 


— 144+0.320 
^^•^^^-0.332 


— 2 fifiq+0.287 


/'(7^o) 


-2.710^16 


1 + 1.8. 10"^^ 
-1.8. 10-15 


1 + 5.8. 10-10 
-5.3. 10-10 


/"(7^o) 


5.9 • 10-20+7.3 10-20 
-7.810-20 


-4.1 • 10-19+4 S. 10-18 
-4. 7. 10-18 


— 1.2 • 10-1®+'^. ^■lo^-'^ 

-6.8. 10-10 


/"'(7^o) 




1.8 ■ 10-9 + 1-610-8 
-1.6. 10-3 


1.8 • IQ-S + i 1 10-8 

-9.5. 10-0 


nno) 






-3.87- IQ-^+f i io";^ 

-8.2-10-12 



Table of numerical references for /(T?.) and its derivatives, evaluated at 2 = 0, i.e. TZ = TZq; the error propagations have been evaluated 

through the standard differential rule. 



Figure 3: 2-dimensional marginalized probability for the pa- 
rameters of model C. The dashing denotes the likelihood of 
the samples. 



ing from data the correct f{Tl) [33] ■ Thus, we get 
1 



fin) 



An' 



2aTTe 



TC/TJo, 



' 2{a + b + c)ennl 
+ e (gB + (a + 2c)7r + 86 arctan 

2no {{a + b + c)7^7^o - 46A) - 



+ 



eR 



{2b - air) An 



2ce7^A(7^ - 7^o)^ sin 



/27^7^^ 

(3i; 



with a, 6, c free parameters of the model. Clearly, with 
this choice for f{n) we obtain /(72.o) = 72-0+ A, /'(T^o) = 
1 and f"{no) = 0, independently of the parameters. 
Next, we calculate the third and fourth derivatives in 
n = 7^0, i-e. 



and 



/ (7eo) = A 



2b + Tr{a- 12c7r) 
{a + b + c)TTnl 



aA 

{a + b + c)nt 



(32) 



(33) 
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Figure 4: 1-dimensional marginalized probability for the pa- 
rameters explored with MCMC. Solid lines (red) are for model 
A, dotted lines (blue) for model B, and dashed (black) for 
model C. The vertical dashed lines are the lower and upper 
limits allowed for the ACDM model by using the WMAP7y 
+ BAO + Ho observations as inferred in [381 ] . Note that for 
the cases /2zo, /szo and f^^o these are very close and cannot 
be distinguised. 



Again, we use equations (jH]) and to write TZq in terms 
of the CS and set the value of A = 2(1— 2(7o)'Ho (according 
to ACDM). Using the numerical values in Tab. IIV[ we 
get the numerical results for the higher order derivatives 
of om model, i.e. 

/"' (7^o) = -5.98 X 10-" 2b + .(a -12c.) ^^^^ 



and 



a + b + c 



3.79 X lO-i^a 



Equation (1311) represents a first example of f{R), satis- 
fying the cosmographic constraints of fCS. We evaluated 
Eq. (|3T|) by using the bounds of Tabs. II and V. We hope 
that such a choice could represent a viable candidate to 
extend the ACDM model limiting case. 

V. FINAL FORECASTS 



In this paper, we addressed the problem of reconstruct- 
ing the correct form of f{Ti), through the use of the so 
called cosmography of the universe. In particular, we 
considered cosmography as a tool to infer cosmological 
bounds on /(z) and its derivatives up to the fourth or- 
der and consequently on f{TZ) and its derivatives, at our 
time. In addition, by considering the class of f{TZ) which 
reduces to ACDM at z ^ 1, we got numerical constraints 
on f{TZ) and its derivatives by relating such quantities to 
the CS. 

Once we rewrite the luminosity distance in terms of the 
filZ) coefficients, we can directly measures them, allevi- 
ating the problems of error propagation. In particular, 
we defined such a set of quantities as the fCS, which can 
be expressed in terms of the well known CS. We found 
the numerical constraints through the use of Monte Carlo 
statistical analyses, by adopting the updated union 2.1 
dataset, the HST bound for Hq and the OHD measure- 
ments. 

In this coarse grained picture, we were able to get strin- 
gent limits for the fCS and we propose a candidate of 
f(TV), able to reproduce the dynamics of the universe in 
accordance with the cosmographic results. We hope that 
the reconstruction of f{TZ) by using the cosmographic 
approach can be extended in future works in order to get 
more relevant constraints on different class of /(72.). 



(35) 
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Appendix A: Luminosity distance in terms of the CS and of the fCS 



In this appendix, we write the formulae for the expansion of the luminosity distance (z) in terms of the CS and 
of the fCS around z = 0. More details can be found in ^2^. The expansion in terms of the CS of the luminosity 
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distance reads 

1 r 



z + 



1 - 90 ) - ^ ( 1 - go + jo - 



+ ^(2 + 5jo - 2go + lOjogo - 15gi^(l + go) + 50)2"* + 



+ 



J__9jo io._Jo_ go_ 11 jo go 27g5 _ 7jogg Hgp _ H^o _ ' 

20 40 12 120 20 12 40 8 8 8 120 

/ J_ Tjo _ 19jo 19^0 , ^ _ 9o 13 Jo go _ 7jggo ^ 7?ogo _ ^ ^ 133jogg 

V30 30 72 720 720 30 9 18 

3 1QQ„4 01„5 1Q„_ Vo_„_ 10„_^_ 7^2 



13g3 



7jogo 133(7f| 21^^; 13so 7joSo , 19goSo , 7ggso 



48 



16 



90 



144 



48 



24 



240 10 



48 



(Al) 



which is a result evaluated at fc = 0; for extensions see [21|. Inverting the system of Eqs. (|25|) to obtain the CS in 
terms of the fCS, we can rewrite Eq. (jAip in terms of the fCS only. We have 



1 r 
Ho 



fo + 2^0 2 
4Hg 



^ , 9/g + 2(36/o - ^o)nl + 108Hg 3 
^ + 72^4 ^ + 



-45/3 + 18/o(-32/o + f,o)nl - 4(567/o - 21/,o + /2.o)H^ - 2592^^ 



"z^ + 



17280-H^ 



576Hg 

945/^ + 2/2(8235/0 - 27Af,o)Hl + 36(2853/^ - 141/o/.o + /'o + 4/o/2.o)H^ + 



24(11151/0 - 459/,o + 3O/2.0 - hzo)nl + 241056H| 



1 



207360-Hf^° 



-8505/0^ + 2/o3(-93555/o + 3214/,o)H^ - 4/o(398115/o' - 22252/o/„ 



462/o/2.o)H^ - 24(271161/2 + /,o(187/,o - IO/2.0) - 3/o(5480/,o - 247/2.0 + 5/3zo))hJ^ 
48(263844/0 - 11478/,o + 843/2.0 - 39/3.o + fizo)n^o - 9315648-^1°) z^l 



(A2) 



